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ABSTRACT 

We study the post-main sequence evolution of products of collisions between main sequence stars 
(blue stragglers), with particular interest paid to the horizontal branch and asymptotic giant branch 
phases. We found that the blue straggler progeny populate the colour-magnitude diagram slightly 
blueward of the red giant branch and between 0.2 and 1 magnitudes brighter than the horizontal 
branch. We also found that the lifetimes of collision products on the horizontal branch is consistent 
with the numbers of so-called "evolved blue straggler stars" (E-BSS) identified by various authors in 
a number of globular clusters, and is almost independent of mass or initial composition profile. The 
observed ratio of the number of E-BSS to blue stragglers points to a main sequence lifetime for blue 
stragglers of approximately 1-2 Gyr on average. 

Subject headings: blue stragglers - stars: horizontal branch - stars: AGB and post-AGB - globular 
clusters: general 



1. INTRODUCTION 

Blue stragglers are main sequence stars that are bluer 
and brighter than the main sequence turnoff in clusters. 
They are expected to have been formed through either 
dynamical processes, such as direct stellar collisions or 
during the disruption of a binary system; or through the 
coalescence of an unperturbed primordial binary in the 
cluster. 

Blue straggler populations have been used to try to 
and cons train the d ynamical evolution of globular clus- 
ters (e.g. lBu3ll993l) . Since the number of collision prod- 
ucts in a cluster should tell us about the rate of collisions 
in the past, it follows that we should be able to constrain 
the density, mass function, or other cluster properties 
from the blue straggler population. In practice this has 
been difficult, since the number of unsettled questions 
about the formation mechanisms of blue stragglers, and 
com peting effects of clust er evolution, is large. How- 
ever, [SmijEii] (|2000D and lMonkman et alj (|2006h claim, 
based on the blue straggler population, that 47 Tucanae 
may have come clos e to core collapse a few Gyr ago; and 
iFerraro et alj (| 19991 ) use the blue straggler population of 
M80 to postulate imminent core collapse in that cluster 
as well. 

But what happens after the blue stragglers evolve off 
the main sequence? Can we see them anywhere in the 
colour- magnitude diagram (CMD), and if we can, does 
that give us any further information about the past of the 
cluster? Blu e stragglers have lifetimes of approximately 
1 Gyr or so (|Sills et al.lll997f ). and it would be useful to 
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probe slightly further back in time if we can. 

There is a population of HB stars in clusters slightly 
brighter than the regul ar horizon t al bra nch that has been 
seen in M3 as early as ISandagel (|1953h . Blue stragglers 
are expected to have masses of up to twice the turnoff 
mass in globular clusters, and using these masses in stan- 
dard single-star evolutionary mod els puts them in this 
part of the HB branch of the CMD (|Renzini fc Fusi Peccil 
1988): BSS progeny should be redder and slightly 
brighter than their low-mass counterparts. Knowing the 
relative numbers of blue stragglers and their evolved 
prog e ny ca n constrain the lifetimes of these two stages. 
Ilbenl (|1986f ) also suggests that different formation mech- 
anisms for blue stragglers could result in very differ- 
ent properties later in their evolution, and so observing 
evolved blue stragglers may help us constrain their for- 
mation. 

Evolved blue stragglers should also exist in the in- 
stability strip, and the expectation is that they could 
be observed as Anomalous Cepheids. These stars are 
low metallicity pulsating stars approximately 2 magni- 
tudes brighter than RR Lyrae stars, and with pulsa- 
tion periods around 1.5 days. Almost all are found in 
dwarf galaxies in the Local Group. They are expected 
to be the result of mass transfer and possibly coales- 
cence in a low mass binary system (up to a bout 1.6 M Q ) 
(|Zinn fc Searldfl976l: ICox fc Proffittj|l98l . which is ex- 
actly the same as one scenario for the formation of blue 
stragglers. Only one anomalous Cepheid has been ob- 
served in a Galactic Globular Cluster: V19 in NGC 5466 
(jZinn fc K ing 1982 ). A detailed spec troscopic study of 
that star ({McCarthy fc Nemed fl997h found a mass of 
1.66 +0.7-0.5 Mq, a small rotational velocity vsmi = 18 
km s _1 , normal iron, s-process, and alpha element abun- 
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dances for a cluster member, and some evidence for a 
very long-period (P ~ 10 000 days) binary orbit. The 
rotational velocity, chemical composition, and orbital pe- 
riod argue against a recent Roche Lobe overflow mass 
transfer episode, but are not inconsistent with a colli- 
sional origin for the blue straggler. Unfortunately, NGC 
5466 is quite a sparse cluster and therefore collisions are 
rather unlikely. Nevertheless, the existence of an unusual 
evolved star as well as a significant number of blue strag- 
glers in this cluster makes it an interesting case. 

The first paper to explore the possible connection be- 
tween blue stragglers an d horizontal branch m orphol- 
ogy in many clusters was iFusi Pecci etall (|1992f ). They 
looked at 21 globular clusters with both blue stragglers 
and well-observed horizontal branches. They argue that 
the reddest horizontal branch stars, particularly in clus- 
ters with bluer HBs, are the descendants of blue strag- 
glers. Based on their analysis of the observations of BSS 
and HB stars in 10 globular clusters, they determine that 
the average ratio of the number of blue stragglers to the 
number of "blue straggler progeny" on the HB to range 
from 1 to 10, with a mean of 6.6 and an overall value of 
5. 3, summed ov er the 10 clusters. 

Bailyn (1994) studied the evolved stars in the nearby 
globular cluster 47 Tucanae, using a combination of 
ground-based and Hubble Space Telescope (WFPC) ob- 
servations to look at the distribution of these stars 
throughout the cluster. He found an enhancement of 
asymptotic giant branch (AGB) stars, approximately one 
magnitude brighter than the horizontal branch, in the 
centre of this cluster only. Since the blue stragglers 
in that cluster (and many others) are centrally concen- 
trated, he mak es the link between th ese two populations. 
More recently, iBeccari et all (|2006fl used HST/ACS to 
look at the evolved stars in 47 Tucanae with better pre- 
cision than the study of iBailynl (|1994l ). They confirm 
Bailyn's result of a central enhancement of AGB stars, 
which they interpret as the result of evolved blue strag- 
glers formed in binary mergers. Neither study was sen- 
sitive to stars much fainter than the HB, and so do not 
give the number of blue stragglers, or the ratio of blue 
stragglers to HB or AGB stars. 

Blue stragglers and their evolved counterparts (E- 
BSS) have been iden tified using HST/WFPC2 in M3 
dFerraro et al.lfl997ah. M13 (jFerraro et al.lll997bl fl999h 



and M80 ( jFerraro et al.|["l999fh The E-BSS are identified 
as being slightly bluer than the RGB and between 0.2 
and 1 magnitudes brighter than the horizontal branch 
in those clusters. In M3 and M80, the radial distribu- 
tion of the blue stragglers and E-BSS is the same, while 
both populations are more centrally concentrated than 
the normal RGB stars in those clusters. The ratios of 
'bright' blue stragglers to E-BSS are determined from 
these studies to be 6.4 in M3, 2.2 in M13, and 6.8 in 
M80 In M80, the population of blue stragglers is so well 
determined that the limitation to bright blue stragglers 
is not necessary. In that case, the ratio of all blue strag- 
glers to E-BSS is - 16. 

In this paper, we use models of stellar collision prod- 
ucts as starting models for stellar evolutionary calcula- 
tions. We follow the evolution of these stars from the 
main sequence to the giant branch, through the helium 
flash onto the horizontal branch, and onto the asymp- 
totic giant branch. We look at their positions in the HR 



diagram and their lifetimes, and compare to the obser- 
vational evidence listed above. We wish to determine 
whether the stars identified in the observational studies 
are indeed evolved blue stragglers, and whether their ob- 
served properties agree with our models of blue straggler 
formation and evolution. We outline our computational 
method in section [21 and present the results in section 
[3j The implications and a discussion of further questions 
follows in section [4] 

2. METHOD 

In this paper, we assume that blue stragglers were 
formed during direct collisions between two unrelated 
main sequence stars. The other methods of form- 
ing blue stragglers (collisions mediated by binary stars, 
and mass transfer in binary systems) will probably 
result in somewhat different structure and composi- 
tion profiles of the newly-bo rn blue stragglers (see e.g. 
BailynjL Pinsonneaultl fl995h . However, the direct colli- 
sion products are the best-studied blue stragglers to dat e 
(|Sills et al.llT997t ISandauist et~aIHl997t ISills et al.|[200l , 
and also the easiest to model. The possible implications 
of this choice on our results will be discussed in section 

SI 

2.1. The collision parents 

We us ed the Yale Rotatio nal Evolutionary Code 
(YREC, iGuenther et all Il992f) to calculate evolution- 
ary tracks for representative globular cluster stars. We 
chose a metallicity of Z=0.001 and Y=0.232 (correspond- 
ing to [Fe/H]=-1.27 where Z Q = 0.0188). Our mixing 
length is a — 1.642, calibrated to give the solar lumi- 
nosity and radius at the solar age for a solar-mass star 
with solar metallicity. We ran models for 0.4 M Q , 0.6 
M Q and 0. 8 Mp, stars, starting a t the deuterium-burning 
birthline (Palla & Stah ler][i991f) . and followed the evolu- 
tion through the pre-main sequence and main sequence 
phases. 

We are interested in collision products which would 
show up in the blue straggler region of the colour- 
magnitude diagram at the current time, or in the recent 
past. We assume that the typical mass of turnoff stars in 
globular clusters is 0.8 M Q . From our evolutionary mod- 
els, the age of a 0.8 M Q star at hydrogen core exhaustion 
is 13.7 Gyr, so we choose this value to be a typical age 
of a globular cluster. 

In previous work (e.g. I Sills et ail I1997D . the parent 
stars were assumed to have the current age of the clus- 
ter. In other words, all collisional blue stragglers were 
formed from parents whose composition profile was that 
of a 13.7-Gyr-old star. However, collision products can 
live on the main sequence for quite a long time - up to 
an additional 10 Gyr for a low-mass collision product of 
0.8 or 0.9 Mq. Therefore, we should be including col- 
lision products whose parents collided quite some time 
ago, and whose structure was that of a much younger 
star. The thermodynamic and physical structure does 
not change much on the main sequence, but the chemi- 
cal composition profile of the star is modified by nuclear 
burning. For the 0.4 Mq models, this effect is not very 
strong since they have completed only a small fraction of 
their main sequence lifetime. However, the 0.8 Mq mod- 
els have significantly different amounts of helium in their 
cores if we look at the models at an age of 10 Gyr com- 
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TABLE 1 

Initial Parameters of Collision Products 



Collision Name 


Mi 


M 2 


Time of Collision 


Mtot 




Mq 


M 


Gyr 


M 


mfMmru A 


0.4 


0.4 


0.03 


0.78 


mll/lmn4R 


4 


4 


2 


0.78 


1 1 1U '-t 1 11 U t: \y 


4 


4 


5 


0.78 


1 1 1W '-t 1 11 L/ t: LJ 


4 


4 


10 


0.78 


1 [ W l<-t- till f~t 1 t 


0.4 


0. i 


13.7 


0.78 


m04m06C 


0.4 


0.6 


5 


0.98 


m04m06D 


0.4 


0.6 


10 


0.98 


m04m06E 


0.4 


0.6 


13.7 


0.98 


m04m08D 


0.4 


0.8 


10 


1.18 


m04m08E 


0.4 


0.8 


13.7 


1.18 


m06m06D 


0.6 


0.6 


10 


1.17 


m06m06E 


0.6 


0.6 


13.7 


1.18 


m06m08D 


0.6 


0.8 


10 


1.37 


m06m08E 


0.6 


0.8 


13.7 


1.38 


m08m08D 


0.8 


0.8 


10 


1.57 


m08m08E 


0.8 


0.8 


13.7 


1.58 



pared to 13.7 Gyr. The younger star has Y=0.855 at its 
centre, compared to Y=0.998 in the older turnoff star. 

We chose 5 different ages of parents to investigate this 
effect: 0.03 Gyr (corresponding to the amount of time 
it takes a 1.0 Mq star to reach the zero-age main se- 
quence from the birthline), 2 Gyr, 5 Gyr, 10 Gyr and 
13.7 Gyr. We labeled these ages as A, B, C, D and E 
when discussing our collision products in the rest of the 
paper. The evolutionary tracks, and the positions of the 
parents at these ages, are show in figure Q] 

In this paper, we looked at all collision products which 
could plausibly be seen in the colour-magnitude diagram 
of a globular cluster - those that could be on the main 
sequence, giant branch, horizontal branch or asymptotic 
giant branch at the current time. To determine which 
collision products to concentrate on, we estimated the 
main sequence lifetime of each product after the colli- 
sion. If the sum of the time of the collision plus the main 
sequence lifetime of the product was 10 Gyr or longer, we 
included those stars in our investigation. For example, 
a collision between two 0.8 Mq stars at time B (2 Gyr 
after the cluster was formed) results in a collision prod- 
uct which is approximately 1.5 to 1.6 Mq. The main 
sequence lifetime of a 1.5 M© star is something like 4 
Gyr, and so this particular collision product will not be 
around at the present time. In fact, the only collisions 
between two 0.8 M Q stars that could possibly be ob- 
served now are those which occurred at times D and E 
in this scheme. Collisions between two 0.4 Mq stars, on 
the other hand, produce stars with masses of 0.8 M Q or 
less, and all of those collision products should still be on 
or close to the main sequence. 

The collisions that we will discuss in this paper are 
given in Table Q] We give the name of each collision, the 
masses of the two parent stars, the age of the parents 
at the time of collision, and the mass of the collision 
product. 

2.2. The collision product and evolution to the Main 
Sequence 

We used the cod e Make Me A Star (MMAS 
lLombardi et al.1 I2002D to calculate the structure and 
chemical profile of our collision products. MMAS takes 
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Fig. 1. — Evolutionary tracks of the parent stars of the collision 
products: 0.4, 0.6 and 0.8 Mq. The stars begin their lives on the 
deuterium burning birthline, contract towards the main sequence, 
and then remain on the main sequence for the full 13.7 Gyr of 
evolution. The properties of the stars at the possible times of 
collisions are marked with large squares, with points A through 
E moving from the pre-main sequence along the main sequence 
evolution. 

the results of stellar evolution calculations (from YREC 
or other codes), and uses a method called "sort by en- 
tropy" to determine the structure and chemical profile of 
the collision product. It has been shown that the results 
of low-velocity stellar collisions can be readily approxi- 
mated by determining the entropy of each shell in the 
parent star (using A = P/p r ) and then creating the col- 
lision product by sorting the shells of both parents stars 
so that entropy increases outward. This technique as- 
sumes that the collision does not shock-heat much of the 
material (thereby changing its entropy), and was based 
on the results of smoothed particle hydrodynamic simu- 
lations of stellar collisions. 

We assumed that all collisions occurred with a perias- 
tron separation of 0.5 times the sum of the radii of the 
two stars. It has been shown that the mass and struc- 
ture of the collision products does not depend strongly on 
the impact parameter, except tha t head-on collisions are 
different than o ff-axis collisions (|Lombardi et al.l 119961 : 
Sills et alJll997h . Since head-on collisions are very rare, 
we chose to use a representative off -axis collision. 

We used the method outlined in ISills et al.l (|1997h to 
take the collisions products and evolve them from the 
end of the collision to the main sequence. For this por- 
tion of the evolution, we used YREC. In this stage, and 
throughout the paper, we neglect any rotation of the col- 
lision product in our evolution. We stopped the evolution 
when the energy generation due to hydrogen burning was 
larger than that due to gravitational contraction, and the 
star had reached the faintest luminosity in that region of 
the CMD. This is equivalent to the definition of the Zero 
Age Main Sequence for a normal pre-main sequence evo- 
lutionary track. This post-collision stage lasted between 
10 7 and 10 s years, depending on the mass of the collision 
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product. 



2.3. Evolution from the Main Sequence 

We wish to follow these collision products on the main 
sequence and giant branch, onto the horizontal branch 
and then to the asymptotic giant branch. These prod- 
ucts are low mass stars, and so they will undergo the 
helium flash at the tip of the giant branch. YREC, like 
most stellar evolution codes, cannot follow the evolution 
of a star through the helium flash because of the very 
short timescales involved and large helium burning lu- 
minosities, up to 10 10 Lq, depending on the initial mass 
and Z. 

The Monash stellar evolution code, howe ver, can evolve 
stars through the core helium flash (e.g., iKarakas et al.1 
2002) 1 . This is because the code has been specifi- 
cally adapted (|Lattanziol 119861 119911 : iFrost fc Lattanziol 
119961 ) to model the evolution of AGB stars and the 
conditions during instabilitie s of the He-burning shell 
(jSchwarzschild & Harm 1965f) are similar, though not as 
extreme, as found in the cores of low-mass stars during 
He ignition. 

Neutrinos losses are also an important component 
when modeling the core helium flash because they move 
the temperature maximum, and hence the ignition point, 
outward from the stellar center. Neutrinos losses in- 
cluded in the M onash code are described in detail in 
Lattanzio (1986). Note that the maximum stellar mass 
that experiences the core helium flash decreases as a func- 
tion of metallicity; this is because these models are hotter 
owing to a lower opacity; at Z = 0.001 this mass is ~ 2 
M Q compared to 2.25 M at Z = 0.02. We stress 
that although we can use the Monash evolution code to 
evolve models through the core helium flash, the mod- 
els should only be considered as crude approximations, 
given that the details of the evolution depend critically 
on the assumptions made about convective energy trans- 
port, that is clearly not one dimensional in nature. A 
thorough understanding about this complicated phase 
will require d etailed multi-dimensio nal hydrodynamical 
modeling, see Dearbor n et al.l (|2006f ) for recent efforts. 

We converted the results of the YREC calculations into 
the format required by the Monash code. We assumed 
that the star was initially fully convective. The Monash 
code determines which shells of the star are convective 
before doing any subsequent mixing or evolution, so this 
assumption was purely for simplicity. To calculate the 
internal energy, we assumed that the equation of state 
was that of an ideal monatomic gas plus radiation pres- 
sure, which is a good approximation to main sequence 
low mass stars. The mass, radius, luminosity, pressure, 
temperature, density and compositions of H, 3 He, He, 
C, N, O and Z were taken from the YREC models di- 
rectly. These new models were then used as starting 
models, and the collision products were evolved from 
the main sequence until the asymptotic giant branch. 
In some cases, the evolution continued through the ther- 
mal pulse phase and onto the post- AGB and white-dwarf 
cooling track. 

In the Monash evolution code we included mass loss 

1 ILattanzTol ||1986D created zero-age horizontal branch models 
instead of evolving through the flash. 
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Fig. 2. — Evolutionary tracks of all the m04m04 stellar collision 
products discussed in this paper, plus the m=0.78 MQnormal star 
track (marked with a dashed line). The other lines are the different 
versions of the same parent stars, collided at different times in the 
parents' lives. 

using the iReimersI (|1975l) formula 

M = -4 x ICT 13 ^ M year" 1 (1) 

where L, R and M are the stellar luminosity, radius and 
mass, respectively, and the parameter rj varies between 
0.4 < n 4, although v alues as high as 10 have been used 
(jStraniero et al.lfl997f l. We set the parameter r\ = 0.4 on 
the RGB and 77 = 1 on the AGB. To test the effect of 
mass loss on the AGB lifetime in the collis i on-pr oduct 
runs, we also used the IVassiliadis fc Woodl (1993) for- 
mulation on the AGB, that relates the fundamental pul- 
sation period of the star to the mass-loss rate, and in- 
cludes a superwind phase once the period increases above 
500 days. We assumed instantaneous mixing, and used 
the same mixing length as used in YREC, a = 1.642. 
For more details of the input physics u sed in the Monas h 
evolution code we refer the reader to iLattanziol fl 986); 
IFrost fe Lattanziol (|1996h and and IKarakas fc Lattanziol 
(|2007ft . 

3. RESULTS 
3.1. Comparison to Normal Stars 

In figures [2] -[71 we present evolutionary tracks of each 
of our collision products. Each figure corresponds to a 
different choice of parent star masses. The different lines 
correspond to the different ages at which the parent stars 
collided. We also present evolutionary tracks of normal 
stars of similar masses to our collision products, shown 
in the figures as dashed lines. The normal stars have a 
uniform chemical composition throughout the star at the 
Zero Age Main Sequence, with Z =0.001 and Y =0.232, 
the same as our parent stars. 

In all these figures, we have marked the approximate 
position of the observed "E-BSS" stars, as well as the 
horizontal branch and the base of the asymptotic giant 
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Fig. 3. — Evolutionary tracks of all the m04m06 stellar collision 
products discussed in this paper, plus the m=0.98 MQnormal star. 
The lines are the same as in figure [2] 
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Fig. 5. — Evolutionary tracks of all the m06m06 stellar collision 
products discussed in this paper plus the m=1.18 MQnormal star. 
The lines are the same as in figure [2] 
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Fig. 4. — Evolutionary tracks of all the m04m08 stellar collision 
products discussed in this paper plus the m=1.18 MQnormal star. 
The lines are the same as in figure [2] 

branch (where these stars are most likely to be found ob- 
servationally) . The HB was defined by the evolutionary 
track of the m=0.78 M Q track, and lies at a luminosity of 
log(L/L Q ) ~ 1.7. The AGB box was chosen to lie just 
above the start of the AGB (as defined below ) for most 
of our tracks. 

The E-BSS stars are, observationally, found somewhere 
between 0. 2 and 1 magn i tudes brighter than the ho rizon- 
tal branch (|Bailvnlll99llFerraro et alJll997aLll999t l. and 
just to the blue of the RGB. Because we are working 
in the theoretical plane, we need to convert from magni- 
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Fig. 6. — Evolutionary tracks of all the m06m08 stellar collision 
products discussed in this paper plus the m=1.38 Mgnormal star. 
The lines are the same as in figure [2] 

tude differences to luminosity differences: one magnitude 
is equal to a difference of 0.4 in logi. Therefore, our E- 
BSS box goes from log(L/L Q ) = 1.78 to 2.1, and we 
chose the temperature extent to begin just blueward of 
the RGB and extend for a few hundred degrees. 

The main sequence position of each track in the CMD 
is very similar in each figure. Choosing a different time 
of collision for the parent stars makes a small difference 
to the position of the collision product in the CMD. The 
evolutionary tracks are also very similar to their nor- 
mal main sequence counterpart. As has been shown in 
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Fig. 7. — Evolutionary tracks of all the m08m08 stellar collision 
products discussed in this paper plus the m=1.58 MQnormal star. 
The lines are the same as in figure [2] 

previous work (|Sills et al.lll997l ). the collisions involving 
lower-mass parents show less effect on the main sequence. 
The lower-mass parents are less evolved overall. The 0.8 
Mq stars are at the turnoff at time 'E', whereas the 0.4 
Mq stars are simply a few more Gyrs along in their ~ 
40 Gyr evolution. 

The post-main sequence evolution of the normal stars 
and the collision products is very similar. There is some 
slight colour difference between the normal stars and col- 
lision products on the giant branch, with the normal stars 
being slightly redder. The difference is usually a few tens 
of degrees at a given luminosity, and is probably not ob- 
servable. The largest different is ~ 300 degrees at the tip 
of the giant branch in figure [5J and even this will be dif- 
ficult to observe since there should be essentially no blue 
stragglers in that evolutionary state in most clusters. 

The location of the core-helium-burning phase for 
higher-mass stars is consistent with the nominal observed 
position of the E-BSS stars in the CMD. For all the 
collision products except the m04m04 series, the E-BSS 
box covers the majority of the track during core-helium- 
burning. 

3.1.1. Horizontal Branch Evolution 

In the CMDs for both normal stars and collision prod- 
ucts, we see some loops in the CMD on the early part 
of the horizontal branch, see figure [5] These loops 
are a product of taking the star through the helium 
flash, rather than stopping the evolution at the tip 
of the giant branch and restarting it with a zero-age 
horizont al branch model, as is done by a number of 
groups (ICharbonnel et all fl9M ICastellani et all Il992t 
IStancliffe et al.ll2005l ). The timescale for these loops is 
very short (< 10 6 years), and they occur at the end of 
the helium flash but before stable core helium burning 
has begun. 

The loops are caused by successive oscillations of the 
core as it settles down after the core He-flash. In fig- 
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Fig. 8. — Evolutionary track of the m04m04B model from the tip 
of the giant branch to the start of quiescent core helium burning. 
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Fig. 9. — (A) Surface (black) and helium-burning (red) luminos- 
ity as a function of scaled time for the m04m04B model from the tip 
of the giant branch to the start of quiescent core helium burning. 
Time is scaled such that the beginning of the flash is t = 0, with 
time in units of 10 6 years. Most of the energy produced by the core 
helium flash (maximum of ~ 10 10 L© at first flash) goes into lift- 
ing the electron degeneracy, with some converted into mechanical 
energy to expand the star. (B) Logarithm of the effective temper- 
ature. The time from the flash peak to quiescent core He-burning 
is approximately 10 6 years. 

ure OJa) we show the surface and helium-burning lumi- 
nosity from the m04m04B collision model during the core 
helium flash. Illustrated are successive flashes of lessen- 
ing intensity, these occur because the core helium flash is 
off center to begin with, and each successive flash mov- 
ing inward (in mass) and removing degeneracy. This is 
similar to the degenerate c arbon ignition found in 8 to 
11 M Q stars (|Siessl I2006D . In figure Mb) we see the 
change to the effective temperature from the tip of the 
GB to the beginning of quiescent core helium burning. 
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The effective temperature first increases when the flash 
begins, followed by short timescale, small amplitude os- 
cillations caused by the whole star expanding and con- 
tracting from the successive core flashes. Once the star 
begins quiescent core helium burning the effective tem- 
perature stays approximately constant. We do not as- 
sume any extra mixing, convective overshoot, etc. from 
the border of the core during this phase. 

3.2. Lifetimes and Possible Evolved Blue Straggler Stars 

We are interested in how many post-main sequence 
blue stragglers we might expect. One easy way to ap- 
proximate this is to look at the timescales of each of the 
different stages in the CMD, and to equate the relative 
number of objects in each stage to the ratio of the life- 
times. First we need to clearly define each stage. We 
took the "zero age" main sequence to be the start of the 
calculation using the Monash code. At this point, the 
stars have come into hydrostatic equilibrium, and are 
burning hydrogen in their cores. We define the end of 
the main sequence (Terminal Age Main Sequence) to be 
the age at which the central helium mass fraction Y has 
reached 0.99. We are also interested in the age that de- 
fines the beginning of the horizontal branch, which we 
take to be the age at which Y drops to 0.97 after the 
TAMS. The end of the horizontal branch phase occurs 
when the central Y drops to zero, and we also take this 
to be the start of the asymptotic giant branch phase. 

Mass loss terminates asymptotic giant branch 
evolution when the e nvelope mass is reduced to 
w O.OlM (|Blockerll2OOll ). or when q = M core /M tota i > 
0.9 with q approaching unity as t he star evolves to th e 
white dwarf (WD) cooling track (|Schoenbernerl [T979D . 
We choose to define a star to have left the AGB 
track when the effective temper ature, T e g, increases 
by A log T G ff = 0.3 on the CMD (|Karakas fc Lattanzid 
2007). From figure ITOl we can approximate the location 
of the AGB to be log Toff = 3.7. Some models evolved to 
the WD cooling track, so to obtain an accurate estimate 
of the AGB lifetime, we needed to remove the time spent 
during the post-AGB (usually short, < 10 5 years) plus on 
the WD cooling track. For other models, convergence dif- 
ficulties ended the calculation before all of the envelope 
mass was lost and the model did not move away from 
the AGB. For all the collision models, however, q > 0.95 
and all final envelope masses are few times 10 _2 M Q or 
less (maximum is O.O34M ). In particular, all Reimer's 
models had q > 0.98 except the m08m08D model, that 
had q = 0.976 and an envelope mass of 0.016M Q . For 
models that had their evolution halted by convergence 
difficulties, the final envelope mass was small enough 
that we can assume that no more thermal pulses would 
occur, and that the final timestep is an excellent approx- 
imation to the final AGB time. This is further justified 
because the thermally-pulsing phase is itself a small frac- 
tion of the total AGB lifetime. For example, the total 
AGB lifetime of the 1.58M , Z = 0.001 normal model 
is 488.6 Myr, much longer than the thermally-pulsing 
AGB lifetime of only 2.1 Myr. The final envelope mass 
was « O.O3M , which would be lost in much less than 
the time taken between successive shell flashes. 

In figure 1 101 we show a few evolutionary tracks, with 
these important points marked with solid symbols. The 
duration of each phase for all our collision products is 
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Fig. 10. — Evolutionary tracks of a few stellar collision prod- 
ucts. From faintest to brightest, they are: m04m04E, m04m06E, 
m06m06E, m06m08E, m08m08E. The zero-age main sequence is 
taken to be the beginning of these tracks. The terminal-age main 
sequence is marked with squares; the zero-age horizontal branch is 
marked with triangles; and the terminal-age horizontal branch is 
marked with circles. The end of the asymptotic giant branch phase 
is marked with stars. 

given in table and for the normal stars in table [31 
In all cases, the mass loss on the AGB was calculated 
using the Reimer's formalism. The first column of each 
table identifies the star or collision product. The next 
three columns give the main sequence, horizontal branch, 
and asymptotic branch lifetime in Gyr. The final two 
columns give the ratio of main sequence to horizontal 
branch lifetime, and the main sequence to asymptotic 
giant branch lifetime. These last two columns are an 
indication of the ratio of the numbers of those kinds of 
stars, under the assumption of a constant formation rate. 

The main sequence lifetimes of collision products are 
shorter th a n normal st ars (as seen in previous studies 
ISills et all (jl997l 120011) ). with the lower mass parents 
showing a smaller effect. This has already been pointed 
out in the evolutionary tracks, where the main sequence 
extent is shorter for collision products with more evolved 
parents. 

For all our stars, the horizontal branch lifetime was 
very uniform - approximately 1 x 10 8 years, regardless 
of parent star masses or age at the time of the collision. 
However, the main sequence lifetimes of these stars can 
change by 3 orders of magnitude, so the predicted num- 
ber of HB stars for every BSS depends strongly on the 
kinds of BSS that are seen in clusters. All our colli- 
sions that resulted from two 0.4 M stars colliding will 
not show up as BSS, so they can be excluded from the 
comparison with data. Their progeny, as well, will be in- 
distinguishable from the normal HB stars in the cluster. 

The asymptotic giant branch lifetime was equally uni- 
form, at about 1.5 x 10 7 years. The normal stars also 
have the same horizontal branch and asymptotic giant 
branch lifetimes as each other and as the collision prod- 
ucts. Post-main sequence evolution of low-mass stars is 



TABLE 2 

Lifetimes and Population Ratios of Evolved 
Collision Products 



Collision MS 


HB 


AGB MS/HB MS/AGB 


Gyr 


Gyr 


Gyr 



846.3 
855.2 
906.5 
889.7 
883.8 
404.8 
343.4 
276.2 

65.6 

39.8 
146.6 
120.5 

56.8 
9.7 

30.6 
2.8 



TABLE 3 

Lifetimes and Population Ratios of Normal 
Stars 



Star 


MS 


HB 


AGB 


MS/HB 


MS/AGB 




Gyr 


Gyr 


Gyr 






0.78 


15.1 


0.097 


0.016 


156.5 


924.5 


0.98 


6.50 


0.091 


0.016 


71.6 


416.7 


1.18 


3.19 


0.092 


0.015 


34.6 


212.7 


1.38 


1.81 


0.094 


0.014 


19.3 


128.5 


1.58 


1.17 


0.095 


0.013 


12.3 


90.0 



extremely robust and largely independent of initial mass, 
structure and composition profiles. 

The observed ratios of E-BSS to BSS stars are usually 
less than ~ 10. None of the normal star models have life- 
time ratios that are in this range, although the highest 
mass star (m=1.58 M Q ) comes closest. We can exclude 
the m04m04 collision models as well, both because their 
main sequence to horizontal branch lifetime ratio is so 
high, and because we could not distinguish either the 
main sequence or horizontal branch stars as blue strag- 
glers observationally. Most observers restrict themselves 
to 'bright' blue stragglers when constructing this ratio, 
which means stars more than ~ 0.6 magnitudes brighter 
than the main sequence turnoff, or a difference of 0.24 
in log(L/L Q ). The m04m06 collision products can fall 
within this range, and so it is not clear if we should be 
considering them in our comparison with the observa- 
tions. If we do, then the average ratio of main sequence 
to horizontal branch lifetimes is 17.7; if we only include 
the m04m08, m06m06, m06m08, and m08m08 collisions, 
our average better matches the observations at 9.3. We 
note that the observational data for M80, which is of 
higher quality than for the other cluste rs, gives a BSS to 
E-BSS ratio of 16 (jFerraro et al.lll999f ). Because of the 
higher quality, blue straggler selection was not an issue, 
and so the quoted ratio is for all blue stragglers, not only 



the bright ones. Therefore, our calculated average of 17.7 
from all our models is in remarkably good agreement. 

From these lifetime ratios, we can also conclude that 
blue stragglers must have short main sequence lifetimes. 
If they behaved like normal low-mass globular cluster 
stars, then we would expect on the order of 100 main se- 
quence stars for every evolved, horizontal branch phase 
star. Since we see ~ 10—20 BSS for every E-BSS, our cal- 
culations point to higher mass collision parents and a late 
time of collision. Blue stragglers should have lifetimes of 
^1 — 2 Gyr to match the observed population ratios in 
M3, M13, and M80, in accordance with our models. 

In figures [TT] and [121 we show the evolutionary tracks 
for the collisions which occurred at times D and E, re- 
spectively. In these figures, the points are equally spaced 
at 10 7 year intervals. By plotting the tracks in this way, 
it is easier to see where the stars would lie in a cluster's 
CMD. The main sequence is very well populated, since 
it is a long-lived phase. The RGB coverage gets sparser 
towards the top, as the stars' evolution speeds up. The 
core-helium burning phase is reasonably long lived, and 
then there are a scattering of points on the AGB, as ex- 
pected. We see very few significant differences between 
these two CMDs. If we assume that blue stragglers are 
formed at a constant rate in clusters (at least over the 
past ^3 Gyr) , then the CMDs of real clusters should look 
something like these diagrams, and we will not be able 
to tell much about the formation rate except by looking 
at the bright end of the main sequence. 

Clearly we have too many points on the main sequence 
- real clusters have a few tens to a few hundred blue 
stragglers per cluster, rather than the thousands of points 
seen in these figures. We are also showing an artificial 
concentrate of stars along the tracks we have chosen, 
when of course blue stragglers could be formed with any 
mass between 0.8 and 1.6 M Q in this scenario. However, 
it is interesting to note that the E-BSS stars all clump in 
the same region, regardless of mass. Therefore, the pres- 
ence of E-BSS stars in a cluster will give us very little 
information about the masses of their progenitors on the 
main sequence, other than the fact that they are slightly 
more massive than the normal turnoff (greater than 1 
M compared to 0.8 M© i n this work). 

The effect of using the IVassiliadis fc Wood! (|1993L 
hereafter VW93) mass loss during the AGB instead of 
Reimer's results in ~ 10 — 20% increases in the AGB life- 
times for some of the lower-mass collision models (e.g., 
m04m04A-D). This difference disappears in the more 
massive models, where the m06m06 and m08m08 mod- 
els have essentially the same AGB lifetime regardless of 
the adopted mass-loss prescription. The small differences 
in AGB lifetimes support our conclusions that the post- 
main sequence evolution of E-BSS is fairly robust, with 
the numbers of stars on the AGB not overly dependent 
on input parameters such as mass loss. 

The final core masses do depend on the choice of mass- 
loss law used during the AGB, with VW93 generally pro- 
ducing larger final core masses compared to the Reimer's 
formula. For example, the m04m04 VW93 models had 
final core masses « O.625M compared to 0.530M Q for 
the Reimers models. The difference was largest in the 
m04m06 set, where the VW93 models resulted in core 
masses of 0.72M Q compared to 0.58M©using Reimers. 
For the m06m08 and m08m08 sets of models the dif- 
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Fig. 11. — Evolutionary tracks for all collision products for col- 
lisions which occurred at time D, 10 Gyr after the cluster was 
formed. The points are equally spaced at 10 7 year intervals. The 
boxes outline the HB, E-BSS and AGB regions of the CMD. 
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Fig. 12. — Same as figure [TT] but for collisions which occurred at 
time E, 13.7 Gyr after the cluster was formed. The boxes outline 
the HB, E-BSS and AGB regions of the CMD. 

ference had essentially disappeared, with both mass-loss 
prescriptions giving similar core masses of O.68M0, with 
the Reimers producing slightly larger final masses. This 
difference can be explained by noting that the VW93 
formula gives low mass-loss rates (~ 10 _7 MQper year) 
until the pulsation period reaches 500 days, after which 
a superwind begins. In the models (e.g., m04m06) with 
the largest core masses, the superwind phase is never 
reached. The core keeps growing until the envelope mass 
is lost, as in all stars, but the timescale for envelope loss 



is much longer without the superwind. The final core 
masses can be compared to WD mass distrib utions for 
clusters (e.g., NG C 6397 lHansen et all 120071). Accord- 
ing to Fig. 25 from lHansen et al.l ( 2007h . core masses as 
high as 0.72M Q for ~ lMQprogenitors are possibly not 
realistic, although we note that the E-BSS will evolve 
somewhat differently to normal single stars. Deep obser- 
vations of WD luminosity functions in globular clusters 
may be able to tell us something about the early blue 
straggler population in clusters. However, disentangling 
evolved blue stragglers from normal white dwarfs will be 
a challenge. 

4. SUMMARY AND DISCUSSION 

We calculated structural models of products of colli- 
sions between main sequence stars, and evolved those 
models along the main sequence, red giant branch, 
through the helium flash, on to the horizontal branch, 
and through the asymptotic giant branch phase. 

We found that the models of evolved blue stragglers 
lie in the correct place in the CMD to account for the 
observed E-BSS stars identified by a number of groups 
over the past few decades. The lifetimes of these stars are 
consistent with the number of observed E-BSS stars com- 
pared to the number of blue straggler stars. We conclude 
that the stars identified as E-BSS by various authors are 
indeed the progeny of blue stragglers. 

The ages of low-mass stars on the horizontal branch 
and asymptotic giant branch are independent of mass, 
collision history or lack thereof, and initial composition 
profiles. The main sequence lifetimes, on the other hand, 
are strongly dependent on these quantities. Therefore, 
the ratio of evolved to main sequence stars strongly con- 
strains the main sequence lifetimes. Since we see approx- 
imately 10 blue stragglers for every E-BSS star, rather 
than the ~ hundred that are predicted from normal star 
models, we conclude that the main sequence lifetimes of 
blue stragglers are quite short, much shorter than their 
normal counterparts. This result has implications for 
the formation mechanisms of blue stragglers. For exam- 
ple, blue stragglers cannot be fully mixed on formation 
(which extends their main sequence lifetimes) , and in fact 
must have very truncated lifetimes. This also has im- 
plications for determinations of formation rates. If blue 
stragglers are a relatively short-lived phase, they must be 
formed at a high rate. The observational identification of 
the E-BSS stars as high-mass, core-helium burning stars 
must be robust in order to make these inferences about 
blue straggler formation. Any observational work that 
can provide direct evidence about the properties of the 
E-BSS stars will be extremely helpful in this question. 
Also, any observations which can identify E-BSS on the 
AGB could potentially provide more constraints on the 
evolutionary status of blue stragglers, based on the mod- 
els presented here. 

Our results are dependent on our determinations of 
the main sequence and horizontal branch ages of all our 
models. While the main sequence ages are very robustly 
determined, the horizontal branch ages could depend on 
our treatment of the helium flash, or of the border of 
the convective core. However, any of these effects must 
be relatively small. First, our HB ages of normal stars 
are consistent with those found in the literature for low- 
metallicity stars. Secondly, the HB ages of both the nor- 
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mal and collisional models are very similar. Any error 
introduced by an incomplete treatment of the detailed 
stellar physics must be both small and constant across 
our models. 

In this paper, we assumed that all blue stragglers were 
formed through direct stellar collisions. The alternative 
formation mechanism, that blue stragglers were formed 
through binary coalescence, may produce stars with dif- 
ferent structures and chemical compositions. It is cer- 
tainly expected that the main sequence evolution of a 
binary merger product will be different than that of a 
collision product. Post-main sequence evolution, how- 
ever, is quite robust to disturbances early on the main 
sequence. Even for our collision products with different 
initial amounts of helium, or different convective struc- 
tures, the post-main sequence evolution is almost identi- 
cal. There are slight differences towards the end of the 
AGB phase and into the thermal pulses, but between 
the turnoff and the tip of the AGB, the tracks, and the 
timescales, are the same. This suggests that unless bi- 
nary mergers result in a star with a significantly different 
structure or composition, the post-main sequence evolu- 
tion will be the same as presented here. The differences 
in structure or composition will need to be quite extreme 
- for example, if a binary merger product loses a lot more 
mass for the same parent stars, significantly increasing 
the total helium content of the star. There are almost no 
models for blue stragglers formed during mass transfer. 
The only paper to address this issue in detail so far is 



the work of lTian et all (|2006f) . who look at mass transfer 
in systems appropriate to forming blue stragglers in the 
open cluster M67. While they do not compare their evo- 
lutionary tracks to normal stars, the tracks of the blue 
stragglers look quite normal after mass transfer ceases. 
Therefore, we expect that blue stragglers formed through 
binary coalescence will have post-main sequence evolu- 
tion that is very similar to that shown in this paper, and 
will also populate the same position in the CMD. 

On the other hand, the main sequence lifetimes of bi- 
nary coalescence blue stragglers are uncertain. The total 
amount of helium in a binary coalescence blue straggler 
should be the same as that formed from a collision of the 
same parent stars. However, it could be distributed dif- 
ferently in the star, and could affect the main sequence 
evolution. If the lifetimes span the same range as the 
collision products in this paper, the number of E-BSS 
stars in the CMD will probably not distinguish between 
the different formation mechanisms, despite our great- 
est hopes. Further models of binary coalescence and the 
subsequent evolution of those blue stragglers will answer 
that question. 
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